-ODOY rOX UBRARt 

- <*a i .aas 



+ 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 




Approved for public release; distribution unlimited. 



T178075 



SECURITY CLASSIFICATION OF THIS PACE (When Dote Entered) 



REPORT DOCUMENTATION PAGE 



nmM u 
WMl POSTGRACiA’ 



t. REPORT NUMICA 



READ INSTRUCTIONS 
ITT SCHOC6EFORE COMPLETING FORM 



2. GOVT ACCESSION NO 



3. RECIPIENT'S CATALOG NUMBER 



4. TITLE (end Subtitle) 



S. TYPE OF REPORT A PERlOO COVEREO 



Operational Availability with 
Spares Support 



7. AUTHOR! A) 



Limited 






#. 



Master' s 'Thesi s ; 
March 1 977 

Performing org. report number 

CONTRACT OR GRANT NUMBERf#) 



Choi , Boo Woong 



9 . 



PERFORMING ORGANIZATION NAME ANO AOORESS 



10 . 



Naval Postgraduate School 
Monterey, CA 93940 



II. CONTROLLING OFFICE NAME ANO AOORESS 



12 . 



Naval Postgraduate School 

Monterey, CA 93940 “ 

[A. MONITORING AGENCY NAME A AOORE SS(ll different from Controlling Office) IS. 



PROGRAM ELEMENT, PROJECT, TASK 
AREA A WORK UNIT NUMBERS 



REPORT DATE 

March 1977 

NUMBER OF PAGES 

53 

security class. (»i mi, report) 



Naval Postgraduate School 
Monterey, CA 93940 



UNCLASSIFIED 

IS*. OECLASSI FI CATION/ DOWNGRADING 
SCHEOULE 



16. DISTRIBUTION STATEMENT (el thle Report) 



Approved for public release; distribution unlimited 



\ DISTRIBUTION STATEMENT (of the ebetrect entered In Block 20, II different from Report) 



8. supplementary notes 



9. KEY WOROS ( Continue on reeeree elde If neceeeery end Identity by block number) 

Coordinated Shipboard Allowance List (COSAL) 
Incremental availability 
Operational availability 



L ABSTRACT (Continue on reeeree elde If neceeeery end Identify by block member) 

Mathematical express i ons are derived for the operational 
availability of a system having a finite number of standby spare 
components. The availability expressions are then used in a 
decision criterion for the selection of items to be stocked to 
support units on extended missions with no resupply capability. 



DD I JAN *73 1473 EDITION OF I NOV 6* IS OBSOLETE 
(Page 1) S/N 0102*014*6601 I 



UN CLAS S I FIED 

SECURITY CLASSIFICATION of THIS PAGE fRI***i Dot* tntor^) 



1 



fuCU^lTV CLASSIFICATION QF THIS P»GEr^«n En<mr+J 



20. Abstract (continued) 

Various algorithms for determining how best to allocate 
constrained resources so as to maximize system operational 
availability are examined. Comparisons with the current Navy 
shipboard allowance list policy reveal that large potential 
savings in resources or better supply effectiveness could result 
from implementation of the procedures of this report. 



DD Form 1473 
. 1 Jan 73 

S/N 0102-014-6601 



UNCLASSIFIED 

SECURITY CLASSIFICATION of this PACErFhifi Dttm Enfrmd) 
2 



Approved for public release; distribution unlimited 



Operational Availability with Limited Spares Support 



by 



Choi , Boo Woong 

Lieutenant Colonel, Republic of Korea Army 
B.S., Republic of Korea Military Academy, 1963 
B.A., Seoul National University, ROK, 1967 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
March. 1977 



6.1 







ABSTRACT 



Mathematical expres s i ons are derived for the operational 
availability of a system having a finite number of standby 
spare components. The availability expressions are then used 
in a decision criterion for the selection of items to be 
stocked to support units on extended missions with no resupply 
capability. Various algorithms for determining how best to 
allocate constrained resources so as to maximize system 
operational availability are examined. Comparisons with the 
current Navy shipboard allowance list policy reveal that 
large potential savings in resources or better supply effec- 
tiveness could result from implementation of the procedures 
of this report. 
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F(t) 



Fixed cost of repair 

Unit spare cost of ith system for time t 
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u 
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R 

X 

T 

TC(k) 

Y 

0 



Constant replacement repair 

Number of failures for time t 

Operational availability of redundant system 
with repair capability at time t, given k spares 

Replacement time random variable 

Constant failure rate 

Time-to-fai 1 ure random variable 

Total cost of system, given k spares 

(u-X ) 
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I. INTRODUCTION 



Most supply systems encountered in the real world are 
multiechelon in nature and all levels contribute to the 
achievement of the support goals of the system. The most 
important support level viewed with the objective of attaining 
a high level of equipment operational availability is the 
last echelon of support, such as the shipboard echelon of 
support in the Naval Supply Systems. This echelon of support 
is responsible for enabling operating units to be self- 
supporting for reasonable periods of time. 

The operational availability of an equipment for a 
certain period of time is generally conceptualized as "the 
proportion of the time an equipment will spend in acceptable 
states." [8] In connection with this, operational availability 
is often used as a measure of effectiveness in the complex 
multiechelon supply systems. The definition of operational 
availability for both components and systems which is most 
often used is "the ratio of the mean time between failure 
(MTBF) to the sum of MTBF, the mean time to repair (MTTR) 
(perhaps with a further partition of MTTR into the mean 
logistics delay time, supply response time and repair time)." 
[10] The steady state operational availability upon which 
the above definition is based is not correct, however, if the 
number of spare parts is finite. 
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Th.e operational availability of an equipment (or a 
system) depends on various factors, such as the reliability 
of its components, the repair and replacement delays, the 
availability of spare parts and supply response times. The 
Naval Supply System influences the latter two variables 
through its resource allocation decisions and policies. The 
system attempts to achieve high availabilities of spare parts 
at the last echelon of support (ship storeroom) through its 
Coordinated Shipboard Allowance List (COSAL). The present 
stockage policy focuses on those components which, histori- 
cally, have experienced sufficient demand to justify spare 
part support. If a component qualifies to be included in the 
COSAL, its depth is determined as that amount sufficient to 
satisfy, with a specified probability (the same probability for 
all components), all demands which might occur in a given 
interval of time. 

The present COSAL procedure has a great deal of intuitive 
and computational appeal. Since it is demand based and 
historical demand data are available, the COSAL is simple to 
determine provided sufficient resources are available. 

However, the current procedure has some shortcomings: 

A. Satisfying demands with a given probability for a 
given interval of time does not guarantee a high level 
of operational availability. 

B. No attempt is made to look at system configurations 
in determining what components to support. 
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C. When resource constraints are active, the use of 
military essentiality codes and other "fixes" must 
often be made to augment COSAL rules to restrict the 
allowance lists. 

Because of these and perhaps other deficiencies of the current 
COSAL policy, improved COSAL determination schemes need to be 
i nvesti gated . 

This paper develops two probabilistic models for 
computing the operational availabilities of systems with 
spare parts and it examines several resource allocation 
algorithms. It offers an alternative procedure for allocating 
spare stock at the shipboard level which could be directly 
related to the operational availabilities of the weapon 
systems being supported, and which would alleviate the short- 
comings of the present COSAL procedure described above. 

Renewal theory and Markovian analysis techniques are employed 
to obtain the results. 

In Chapter II, the thesis develops a mathematical 
expression for the operational availability of a redundant 
system where there is no repair capability. Mathematical 
expressions for computing increment availabilities by adding 
an extra spare part are also determined in Chapter II. In 
Chapter II another mathematical model is developed for the 
determination of the operational availability of redundant 
systems where a repair capability exists. For modeling 
simplicity throughout the thesis, the mean time to failure. 
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tile mean time to repair, and the mean time to replacement 
are all assumed to have exponential distributions. The first 
section of Chapter IV addresses a decision rule which 
considers the cost of spares and the cost of repairs to 
determine whether an extra spare should be carried or if a 
repair capability should be bought. Another decision rule 
which considers the "cost of system beinq down" as well as 
the cost of spares and cost of repair is examined in the 
second section of Chapter IV to determine the optimum number 
of spares which should be carried. In the last section of 
Chapter IV an algorithm to determine how best to allocate a 
fixed budget is presented. That algorithm ignores system 
configuration, treating each component as a system. The 
algorithm attempts to maximize operational availability per 
dollar expended on spares support. Chapter V shows how the 
resource allocation scheme can be affected by considering 
the system configurations. Each chapter includes numerical 
examples to illustrate the mathematical results. Summary 
and conclusions of the paper are presented in Chapter V, 
together with suggested extensions of this study. 
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II. OPERATIONAL AVAILABILITY FOR A SYSTEM WITH STANDBY 
SPARES AND NO REPAIR CAPABILITY (TYPE-1 SYSTEMj 



Mathematical expressions for computing the operational 
availability of a system having a finite number of spare 
components are derived in this chapter. It is assumed that no 
repair capability exists for the system, so that, after all 
spares are used, the system must remain down if a failure 
occurs. Mathematical expressions are also derived for the 
incremental change in operational availability for additions 
to the number of spares carried. These expressions form the 
basis for the work presented in later chapters. 

There have been many studies and papers written on 
standby redundancy in systems with repair maintenance. See, 
for example, [5, 9, 11, 12]. However, for the problem 
considered here, the systems often cannot be repaired at sea 
because of lack of a repair capability or other resources or 
because of constraints of space and time. Apparently, little 
work has been done for the case of finite spares and no repair 
capabi 1 i ty . 

A. MATHEMATICAL EXPRESSIONS FOR OPERATIONAL AVAILABILITY 

(TYPE-1 SYSTEM) 

Consider a system of a single component having exponen- 
tial times to failure and exponential times to replace with 
mean time between failures (MTBF) and mean time to change 

(MTTC) equal to 1/A and 1/y, respectively. Suppose that k 
identical spares are provisioned and available for support 
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at the start of a mission. When an operating unit fails, 

replacement begins i mmedi ately provided a spare unit is 

t h 

available. Assume that the time to failure of the i unit, 

, and the time to change, are independent, and the 
system under consideration is operational ("up") at time 
t = 0. Finally, suppose that non-operating idle spare parts 
are not subject to failure or deterioration. A system with 
these properties will be called a "Type-1 System." 

For k = °° , the up and down availability are repeated 
over the duration of the total mission of the ship. The 
availability status of the system can be described as an 
alternating renewal process shown in Figure 1. 

| x O x O x — O 1 

Start Down Up Down Up Down Up Terminate 

Mission Mission 

OPERATIONAL STATUS OF SYSTEM WITH UNLIMITED SPARES 

Figure 1 

For this case, the operational availability of the system at 
time t is easily shown [1] to be: 



AJO = 



MTBF 



MTTC 



MTBF+MTTC MTBF+MTTC 



exp ( - t ( 



1 



1 



MTBF MTTC 



)) 



( 1 ) 



On taking the limit as t the steady state operational 
availability is found to be: 



1 i m 



Ajt) 



MTBF 

MTBF+MTTC 



( 2 ) 
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Equation (2) is widely used as a definition of operational 
availability. However, as stated above, it is exact only 
for the case where the number of spares is unlimited. If 
the number of spare parts is finite, the operational status 
of the equipment is not described by an alternating renewal 
process. Once the spares are exhausted and the system goes 
down, it will remain down for the remaining duration of the 
ship's mission or perhaps the ship's mission must be aborted. 

Of course, if the number of spares is large or if the length - 



good approxi mati on to the true operational availability. In 
times of tight resources such as those currently experienced, 
it is not possible to maintain large pools of spare parts. 
Therefore, it is important to know just what is the avail- 
ability for any finite number of spares. Resource allocation 
algorithms must know the benefits to be derived from 
incremental spare part allocations in order to allocate 
resources efficiently. Thus, expressions for the availability 
for any finite number of spares must be determined. 

Define the following probabilities: 




A^(t) = P^Csystem is up at t | k spares and system 
is up at time 0] 



r 



B k (t) = P r [system is up at t | k spares and system 
is down at time 0] o _ 



VC 




t 




V- 
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Up 



Down 



0 



x 



t 



T i me 



Start 
Mission 
(k spares) 

Figure 2. Operational Status with k Spares (Up at time 0) 

Figure 2 illustrates a single cycle of the system with a 
finite number of spares. Now, it is important to consider 
as the "state" of the system not only its operational status, 
but also the number of spares available. The general formula 
for the operational availability at time t can be obtained by 
conditioning on the lifetime of the component, T^. This 
yields the expression: 



where F(x) = P [T<x], (The system is up at time t if the 
first unit survives to time t or if the system goes down at 
time x<t, and the system is back up and operating at time 
t.) Expression (1) can also be written as 




o 



(3) 




o 
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where H(x) = P r [T + R<_x] and M^(x) is the renewal function of 
H(x). Now, let us see how B ^ _ i ( t ) can be expressed in terms 
of A| < (t). Suppose the system is down at time 0 and is 
restored to up condition at time x. This is shown in 
Figure 3. 



Up 



Down 



t-x 



Start 
Mission 
(k spares) 



T i me 



Figure 3. Operational Status with k Spares (System 
Down at Time 0) 



Then, conditioned on R-j = x. 



( 0 for x>_t 

B k (t)/R 1 = x = ) 

( A k _ i ( t-x ) for x<t . 



Therefore , 



B k (t) 




A k-1 (t-x)dG(x) 



(4) 



where G(x) = Pr[R<_x]. Substituting (4) into (3) gives 



A k (t) 



F(t) + 




X 

A k _ 1 (t-x)dG(x)]dF(x) 



(5) 
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This can be simplified by taking the Laplace transform of 
both sides giving: 



A k (s) = F(s) + B k (s) • F ( s ) 



and from (4) B k (s) = A k _-|(s) * B ( s ) 



Substituting (7) into (6) gives 



A k (s) = F(s) + A fc _ (s) • H ( s ) 



Inverting (8) gives 



A k ( t ) = F(t) + I A k _ 1 (t-x)dH(x) = 



r 



F(t) + J A k _i (t-x)h (x)dx 



( 6 ) 

(7) 

( 8 ) 

(9) 



where h(x) is the density of H(x). Equation (9) is a general 
formula for computing the operational availability of a 
Type-1 System at an arbitrary time t for k = 1,2,*-', 

For k = 0, it is clear that A Q ( t ) = F(t) = e" At . The expressions 
obtained from (9) for k = 0,1,3 and 4, and t >_0 , are summarized 
below with 0 = Xy and y = y-x: 



A Q ( t ) = e' Xt 



( 10 ) 



A i ( t ) = (1 



0 , 0tx -Xt 

~~0 + ) e 

2 y 

Y 



+ 



e -yt 

2 e 
Y 



(ID 



A 2 ( t ) = [1 






Y 




]e 



-xt 




3q 

2 

Y 

1 7 



* 




- vt 




( 12 ) 



(13) 
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(14) 



The derivations of the exact expressions from equation 
(9) get tedious rapidly as k gets large. Consequently, an 
alternative approach was tried. Write (8) recursively for 
k = 1 , 2 ,..., n to obtain: 

A^Ts) = yis) + • 7(T) = F ( s ) [ 1 +700] (15-1) 

A 2 (s) = R7) [1 +7U) + Hit) 2 ] (15-2) 

v"s) =?(t) [1 +'!?(t) + hTT) 2 ( 1 5-n ) 

+ H (s ) 3 +•••+ H (s ) n ] . 
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On inverting and simplifying algebraically we can determine 
the solution. For example, when n+«, ( 1 5 - n ) becomes: 



limA (s) = F(s) [- 



1 



n+oo 



l-H(s) 



-] 



(16) 



since H(s) < 1 and a geometric series exists. Thus 

V 

~ J S 1 n _ / S x/ 1 



Tim A n (s ) = (—rr ) C 

n^oo 



= i - 



XS 1 -F ( s ) G ( s ) 
X 



] = Ht>< 



i 



Xu 



U + s) (y+s) 



s+ ( X + y ) 



(16a) 



Inverting gives 



J r\ 



1 i m A n ( t ) = | 6(x)dx-x 

J o v o 



g- ( U + X ) 



_ H + A - ( u+x ) t 



y + X y + X 



(17) 



where 6(t) is the unit impulse function [4]. Taking the 
limit as t-*-« gives the steady state operational availability 
identical to Eq. (2) as it should be. Unfortunately, the 
inverse Laplace transforms of the expressions for A k (s) 
are not in general easily obtained. 



B. INCREMENTAL AVAILABILITY AND NORMAL APPROXIMATION 

Now, we seek to determine expressions for the incre- 
mental benefit that can be obtained by adding an additional 
spare. Let A ^ ( t ) be the incremental system availability 
that will result from the addition of one spare unit to a 
system already having k-1 spares. Then, 
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(18) 



A k U) = A k (t) - A k-] (t) 

For small values of k, (18) can be evaluated from (9) by 
computing A k (t) and A k _^(t) and s ubstracti ng . When k gets 
large (3 or more) the computations become tedious, as noted 
earlier. Therefore, some approximation method is needed to 
estimate (18) for large k. Using (15-n), we can write (18) 
in terms of its Laplace transforms as: 

Qt) = F(s) • H ( s ) ( 19 ) 



and, inverting. 



A k (t) 




[1 - F(t-x)]dH k (x) 



- H k (t) 




F(t-x)dH k (x) 



( 20 ) 



However, note that H ^ C t ) = P [T.|+R.j + T 2 + R 2 + *** + 

and 
t 

F(t-x)dH k (x) = P r [T 1 +R ] + VR 2 +---+T k + R k + T k + l- t] * 

For large k, we can appeal to the Central Limit Theorem to 
use the normal approximation for the distribution of the sum 
of 2k random variables. Thus, 

H k (t) « N[k(} + I), kt-’-j + ij-)] 

X U 
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and 



I F(t-x)dH k (x)~ N[(k(l + I) + 1). (kflj + Ij) + iy)]. 

J r\ 



Hence, for large k, (18) can be approximated by 



A k (t)~Aj;(t) = *[ 



t_k (T + T) 

i-L] - ♦[ 






] 



( 21 ) 



where 



CT k 




and a k 1 




1 
v 



) + 




Table 1 gives numerical expressions for A k (t), A k (t) and 
A k ( t ) for various values of k and the parameters MTBF = 30 
days and MTTC = 5 days. 

Table 1. Sample Calculations for A k (t), A^(t) and A k (t) 
(MTBF = 30 days, MTTC = 5 days) 





t 


45 Days 






90 Days 




k 


A k (t) 


A k (t) 


*k<‘> 


A k< t > 


A k U) 


A^(t) 


0 


0.2231 


N/A 


N/A 


0.0498 


N/A 


N/A 


1 


0.5712 


0.3481 


N/A 


0.2171 


0.1673 


N/A 


2 


0.7784 


0.2072 


0.133 


0.4622 


0.2451 


0.254 


3 


0.8429 


0.0646 


0.064 


0.6721 


0.2098 


0.159 


4 


0.8554 


0.0125 


0.027 


0.7909 


0.1188 


0.086 


5 






0.012 






0.045 


6 






0.006 






0.023 


7 






0.003 






0.012 


00 








0.8570 


0.000 


0.000 



(1) N/A: Not Applicable (2) A^(t) is normal approximation 
for A k (t) . 
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Figure 4 shows a plot of A ^ ( t ) for k = 0,1, 2, 3 and 4 for the 
parameters MTBF = 30 days and MTTC = 5 days. The curves 
show that little incremental benefit can be gained by 
increasing the number of spares beyond k=4. 

The approach described here would appear to be fruitful. 
For small k, the expressions could be evaluated exactly and 
for k>4, the normal approximation could be utilized. 
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Figure 4. Operational Availability, A^(t). 
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III. OPERATIONAL AVAILABILITY FOR A SYSTEM WITH STANDBY 
SPARES AND A REPAIR CAPABILITY (TYPE-2 SYSTEM) 



In addition to providing a larger number of standby 
spares, the system availability can be imoroved for a given 
spares allotment by increasing the equipment reliability or 
decreasing its replacement time or adding an onboard repair 
capability. This chapter investigates the latter alternative. 
The objective is to derive an expression for computing the 
operational availability of a system with standby spares 
with a repair capability. The tradeoff problem of adding 
more spares or adding a repair capability will be investigated. 

Consider a system consisting of a single component that 
is subject to failure. The system is supported by a finite 
number k of standby spares which, when idle, do not 
deteriorate. When the system fails, a standby spare is 
instantaneously switched on (if a spare remains) and repair 
of the failed component commences immediately. (This implies 
that the replacement time is negligible.) We also assume 
that the repair facility has ample servers so that no failed 
unit has to wait for repairs to begin. (The repair facility 
is a (k+1) server queue.) Under these assumptions, the 
system (called a Type-2 System) is down only when the system 
fails and no standby units are available. Assume that the 
system fails and no standby units are available. Assume that 
the system and all standby spares are up at the start of the 
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mission. As before, assume that the times to failure for 
units in operation are exponential with parameter A, and the 
times to repair are exponential with parameter 3 . Finally, 
assume that the times to failure and the repair times are 
i ndependent . 

Let P.j(t) be probability that i units are in repair at 
time t for i = 1,2,..., k+1. Then, the operational avail- 
ability of a Type-2 System with k spares is 

Q k (t) = 1 - P k+1 (t). (22) 

For the special case where k=0, the system goes down 
and remains down for a repair interval every time the system 
fails. As soon as the repair is completed, the system 
becomes operational. Therefore, the process acts as an 
alternating renewal process, just as with the Type-1 System 
with k=». By analogy, the operational availability is given 
by : 

«o (t) * XT? + IF? 6 , Pt- < 23a > 

and the steady-state operational availability is 

lin q 0 (t) = ^ . ( 23 b) 

t-^oo 

For the general case of k>0 spares, a different 
approach is taken. Assume first that k=l. Then, in a 
small interval of time At, the transitions described in 
Figure 4 can be made. 
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P 2 ( t+A ) 

P 1 ( t + A t ) 
PgU + At) 



Figure 5. Transition Diagram (Type-2 System). 



1 = P (no repair occurs)* e~^ A ^ ^*l-2gAt [7] 

2 = P r (1 unit gets repaired) ^ 2gAt 

2 2 

3 = P r (2 units get repaired) ~ g (At) ^ 0 

4 = P^ (active unit fails and no units are repai red ) ss XAt 

5 = P^ (active unit stays up and no units are repaired) 

+ P^ (active unit fails and a failed unit gets 
repaired) = l-(x+g)At 

6 = P r (active unit stays up and a failed unit gets 

repai red ) gAt 

7 = P (simultaneous failure of 2 units )c^0 

r < 

8 = P^ (active unit fails) ~ XAt 

9 = P r (active unit stays up) srl-XAt. 

Then, the following 3 simultaneous equations hold: 

! P 0 ( t+A t ) = ( 1 - X A t ) P 0 (t) + g A t P -j ( t ) + 0 . P 2 (t) (24) 

P 1 ( t + A t ) = XAt P g ( t ) + [ 1 - ( X + g ) At] P-j(t) + 2gAtP 2 (t) 

P 2 ( t + A t ) = 0 * P q ( t ) + XAtP-j(t) + ( 1 -2gAt ) P 2 ( t ) 
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Rearranging (24) and taking the limit as At+«, yields the 
3 differential equations: 



1 i m 
A t->-0 

1 i m 
A t-*-0 



P 0 (t+At)-P 0 (t) dP Q (t) 



1 i m 
A++0 



At 


dt 


P i ( t + A t ) - P^t) 


d P -j ( t ) 


At 


dt 




+ 2 6 P 2 ( t ) 


P 2 (t+At) - P 2 (t) 


dP 2 (t) 


At 


dt 


solve ( 25 ) , we 


take the L 



= - XPq ( t ) + 6 P i ( t ) 



= xP 0 (t) - ( x+e ) P^t) 



AP 1 Ct ) - 2 6 P 2 ( t ) 



( 25 ) 



sides of each equation and convert the transformed system of 
equations into a matrix representation. The above 2 steps 
are shown in (26) and (27) respectively. 



s P Q (s) - 1 - - X P 0 ( s ) + &P i ( s ) 



s P 1 Cs ) 



s P 2 (s) 
dP,(t) 



PqCs ) - Cx+e)P 1 (S ) + 2 bp 2 (s) 



XP 1 (s ) - 2BP 2 (s) 



(26) 



Since — ] = s P^(s) - f(0~) holds [6], and the initial 

value of P.j(t) is equal to 1 in this problem (say, 

P 0 (0) = f(0“) = 1), (26) holds. 
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( 27 ) 



s+x -8 0 




-P 0 (s)- 




_ 1 “ 


-X s + ( X + 8 ) -28 




p 7 <s) 


= 


0 


0 -X s + 2 8 




l 

to 

1 




0 



The solution for P 2 ( t ) is given by 

V‘> - SF n + TITtrr C b e' at - ae- bt ]}, ( 28 ) 

where a = ( 3^+2 \ ) -V B 2 +4> 6 

_ 3B+2A +V S 2 +4XS 
‘ 2 

Because P 2 (t) is the probability that the system will be 
"down" at time t for the case where only one spare is 
available, the operational availability of a Type-2 System 
for k = 1 is 

Q 1 ( t) = 1-P 2 (t) 

■ SF {1 + T^FT Cbe- at -ae- bt ]}. (29) 

We are also interested in the asymtotic probability of (29), 
which is given by 



11m Q 1 (t ) 



2e(x+e) 

2e 2 +2xe+x 2 



(30) 
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Analogous to the simple case above, for an arbitrary k, 
we have a (k+2) x (k+2) matrix system for solving P k+1 (s): 



(s+x) 

-x 

o 

o 

0 

1 
« 
i 
i 



-e 

(s+x+3) 

-x 

o 

o 

i 

i 

i 

i 

o 

0 



0 

-23 

( s+X+26 ) 

-X 

0 

i 

i 

i 

i 

0 

0 



0 - . 
0 - 
-36 

(s+X+33) 

-X 



0 

0 



-43 - - 

( s + X+46 ) - 



- 0 

- 0 

- 0 

- 0 

- o 

I I I I I 

I I I I I 

I I III 

I I I I I 

-X (s+x+ks ) -(k+1 )g 

- - - -X [ s +x + ( k I 



•1 






“ * 




P 0 (s) 




1 




pf(s) 




0 




P^(S) 


= 


0 




£ (s) 




0 




P 4 (s ) 




0 








♦ 

m 




• 




♦ 




V s > 




0 


] 


>,<•> 




0 

L'J 



Solving ( 31 ) 



for P k + 1 (s ) gives 



(k + 1) 



k+l ^ s[ (s + a 1 ) (s + a 2 ) - - - (s + a k+1 )J 



(32) 



where -a-| > -a 2 j ~ ~ ~ * _a k + l are zeros °f tlie polynomial 

in the denominator. All of the roots are functions of X and 
B , and they can be found numerically by using various computer 
programs. When the roots are known, (32) can be expanded 
using partial fractions and then converted to find P k+ -|(t). 
This gives 



P 



k+1 



(s) 




s s + a.| s + a 2 



k+1 



s + a 



k+1 



(32a) 



( 31 ) 



28 



are all constants and 



where d 1 , d 2 , - - - , d 



k+1 



™ a 1 1 ■ d p t -a.t a k + 1 1 

P k + 1^ = d l e ” d 2 e " d 3 e ' d k + l e 



Then, the general formula for computing the operational 
availability for a Type-2 System is 

" d -I t - d ~ t 

Q k ( t ) = l~(d-|e ~ ^ 2 e “ ^3 e - - - + ^ e 



a k+l t ) 

(33) 



For example, when k=2, the operational availability function 
is given by 



Q 2 (0 1-x ^ ai a 2 a 3 ' a~j (a 2 -a 1 ) (a 3 -a 1 ) ex P( _a i t ) 



(34) 



- a 2 (a r r 2 T(a 3 -qy exp( - a 2 t) - a 3 (a r a 3 )Ca 2 -a 3 )- exp( - a 3 t)] 



Suppose that MTBF = 30 days, and MTTR - 5 days. Then, 

p _ 0.000037 

r 3 ls; ‘ s ( s + 0. 20033 ) (s + 0. 41 09 7 ) (s + 0. 68870) ’ 



and 



P 3 ( t ) “ 0-000037 C( 0 . 20033) (o! 4 1097 ) (0. 68870 ) ‘ 

-0. 20033t 

6 

' TO 0033 ) (0.4T097 -0.20033) (0.68870-0.20033) 

-0. 41 097t 

‘ (0.41097) (0.20033-0.41097) ( 0.68870-0. 4 I 09/)" 

e -0.68870t 

" (0.68870) (0.200 33-0.68870 ) (0.4 1 097 -0.o88'7T3)^ 
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and 



Q 2 ( 45 days) = 1 - P 3 ( 45 days) = 0.99935. 

Various computed values of Q k ( t ) for k = 0 through k = 4 are 
tabulated in Table 2. 

Table 2. Tabulated Values of Operational Availability 
for Type-2 System (MTBF = 30 days and 
MTTR = 5 days) 



X 


45 days 


90 days 


00 


1 


0.85715 


0.85714 


0.85714 


2 


0.98824 


0.98823 


0.98823 


3 


0.99935 


0.99935 


0.99935 


4 


0.99998 


0.99998 


0.99998 


5 


1.0 


1 .0 


1 .0 
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IV. REPAIR/SPARES TRADEOFF AND RESOURCE ALLOCATION ALGORITHMS 



In the last two chapters we have shown that a system's 
operational availability can be improved by either increasing 
the number of spares dedicated to the support of the system, 
or by purchasing an onboard repair capability for the system. 
Both of these solutions involve difficult resource allocation 
decisions. In the case of adding more spares there is the 
problem of finding sufficient space to carry the support units. 
Also, because of tight procurement budgets, the various 
systems must compete for the funds. This creates a real 
resource allocation problem. 

In the case of purchasing a repair capability there are 
perhaps even more difficult problems of manpower planning and 
training; purchasing and maintaining test equipment and repair 
equipment, and, also, the purchase and storage of spare stock. 
With the reduced manning levels that are being imposed on 
ship crews and with the sophistication of weapon systems, the 
recent trend has been away from onboard repair and more toward 
the removable module/replacement alternative. However, a 
tradeoff does exist and it may be that some systems should be 
Type-1 and others Type-2. In this chapter, we cannot quantify 
all of the costs associated with purchasing a repair capability, 
nor can we consider all of the actual resource constraints. 
Instead, we assume a simple cost structure for each type of 
system and look at the cost-benefit tradeoff decision. 
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Later in this chapter we consider the resource alloca- 
tion problem. Given a repair/no repair decision has been 
made about each system and a fixed procurement budget, how 
should the dollars be spent among the many systems that 
require support? We also show the dramatic effect that the 
system configuration can have on the allocation procedure 
and on system availability. 

A. DECISION RULE FOR REPAIR/SPARES TRADEOFF 

Let us assume that the system consists of a single 
component that is subject to failure. Assume that the total 
expected repair cost for a mission of length t is linear, 
i . e . 

C f +C r E[N(t)] 

where N(t) is the random number of failures during the 
mission of length t, C r is a cost (average) per unit of 
repair, and is a fixed cost of repair that is independent 
of the number of units actually repaired. would include 

the cost of training repair personnel, the cost of manning 
for repair, the cost of test equipment, etc. Further, 
assume that an operational availability goal for the mission 
is given. Then, our objective is to decide whether or not 
the repair capability should be purchased and what number of 
spares should be carried. The criterion is the minimization 
of costs. Stated as a nonlinear program: 
jmin(cost) 

( s .t. A k (t)[oriQ k (t)] >_ p 
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The spare stock procurement cost is simply C p -k where C p is 
the procurement cost per unit and k is the number of units 
stocked. The solution to (35) can be obtained as follows: 

(1) Let k be the minimum k such that A k ( t ) >_ p. 

★ k 

(2) Let k be the minimum k such that Q ^ ( t ) >_ p. 

(3) Compare C p *k with C p *k + (C^+C r E[N(t)] } 

t 

where N(t)^rx J Q k **(t)dt 
o 

k k k 

(4) a. If C p • k > C p • k + {C f +C r E[N ( t ) ] } , choose the 

★ ★ 

Type-2 system with k = k 

•*> -k k 

b. If C *k <C *k + { C .p+C E[N ( t ) ] } , the choose 

p p t r 

★ 

the Type-1 system with k = k . 

c. If the costs are equal, select the maximum of 

Q k** ( t ) and 

This decision problem is demonstrated with the following 
numerical example. 

Suppose that p = 0.85, t = 90 days, MTBF = 50 days, 

MTTC = 1 day, MTTR = 5 days, C p = $400, C f = $500, and 

C y = $200 per repair. 

Step 1 . A 3 (90) = 0.8859 -> k* = 3 

* ★ 

Step 2 . Q q (90) = 0.9090 -*■ k =0 
Step 3 . C p *k = $1,500 

E[N(t)]~x f Cxfi- + V5T e- (x+8 > x ] dx 

E[N(90)]s=l .6363 

k k 

C p k + (C f +C V E[N(90)]} = $727.4 
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Step 4 , C p k* >C p k** + (C f +C y E[N (90)] } 

Type-2 System with k=0 is preferable. 

B. DECISION RULE FOR THE OPTIMUM NUMBER OF SPARES FOR SUPPORT 

The problem of realistically allocating spares for the 
support of equipment should take into consideration the "cost" 
of the system being down as well as the cost of spares and 
the cost of repairs. However, the problem of quantifying 
the "cost" of a system being down is very difficult. The 
difficulty lies in the fact that the cost is not so much in 
monetary units, but more in terms of reduced effectiveness 
where the word effectiveness is ill defined. Military 
planners have long agonized over the problems of measuring 
effectiveness and quantifying the cost of a system being down. 
We are still far from a solution today and the solution is 
beyond the scope of this paper. Nevertheless, as in the 
previous section, we assume a simple cost structure with 
known costs for the purpose of demonstrating a solution 
methodology. 

Assume that the cost of the system being down is 
proportional to the down time and that the cost rate is 
known. Also, assume that the total expected system cost 
(TC) is the sum of the cost of spares, the expected cost of 
repairs, and the expected cost of the system down time. 

(For the Type-1 System the expected repair cost is omitted.) 

The cost of spares and the cost of repairs are as defined 
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in the previous section. Thus, the total expected cost of 
system downtime for the mission of length t is: 

C d -E[D(t)] 

where D(t) is the down time in (0,t), and the total expected 
system cost with k spares is: 

TC(k) = C p ♦ k + {C f +C r -E[N(t)]> + C d * E[D (t ) 3 • 

for a Type-2 system and 

TC(k) = C p -k + C d -E[D(t)] 

for a Type-1 system. 

With this cost structure, the algorithm for determining 
the optimal number of spares for either of the two types of 
systems is: 

( 1 ) Let k = 0 

(2) If TC(k+l) - TC(k) >_ 0 set k =k and stop. 

(3) Set k=k+l and go to (2). 

Without additional examination of the costs as functions of 

k, no general claim for optimality can be made. However, 

actual experience has shown that such simple algorithms 

usually yield optimal or near optimal solutions. A numerical 

example of the procedure is given below. Suppose that 

C d = $100/day and that all the other costs and parameters 

are identical to those of the example in the previous section. 

★ 

We determine the number of spares, k , for both Type-1 and 
Type-2 Systems. 
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(a ) Type-1 System : 

TC ( 1 ) - TC ( 0 ) = - 2155.79 
TC ( 2 ) - TC91 ) = - 793.49 

TC (3 ) - TC ( 2 ) = 3 . 31 > 0 - k* = 2 

(b) Type-2 System : 

TC ( 1 ) - TC(O) = - 249.83 

TC (2 ) - TC ( 1 ) = 465.23 > 0 -*■ k* = 1 

C. RESOURCE ALLOCATION SCHEME FOR INDIVIDUAL SYSTEM 

Now, we consider the resource allocation problem. The 
budget available for supporting all the spares that are 
needed will not usually be sufficient to achieve the levels 
of operational availability that are specified for the 
systems. The budget constraints force us to consider 
effective allocation schemes through which a fixed number of 
dollars (B) can be allocated to the many systems being 
considered so as to maximize the operational availability per 
dollar spent on spares. This can be done by marginal analysis 
using benef i t-to-cost tradeoffs. In this section, an algo- 
rithm for the selection of spares for support is developed. 

The algorithm considers the marginal increase in system 
operational availability per dollar spent. With such an 
algorithm a list of spares could be determined which repre- 
sents the most economical way that supply can achieve its 
availability goal. 
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Suppose there are n independent Type-1 Systems. Define 
the following expressions: 

(a) P(-j) = the specified operational availability 

goal for the ith system. 

(b) C(-j) = the unit cost per spare associated with 

the ith system. 

(c) A^(t) = A k (t) for the ith system 

(d) A ^ 1 ^ ( t ) = A k ( t ) for the ith system. 

Then, the recommended allocation algorithm is as follows: 

Step 1 . Compute A^ ^(t) for i = 1,2,3, ... , n. 

If A^(t) >_ for every i, then no 

allocation of dollars is needed. Stop. 
Otherwise, go to Step 2. 

Step 2 . If B < Min {C^} for ith system such that 

A ^ 1 ^ C t ) < P(-j)> then allocation is impossible. 
Stop . 

If B = Min {C^j} for ith system such that 
A^(t) < P^), then allocate C^. Stop. 

If B > Min for ith system such that 

A<J>(t) < P(-i)’ then go to Step 3. 
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Step 3. 



Compute Max { 



} (the greatest increase 



C (i) 

in system availability per dollar spent). 

Suppose the maximum occurs when i=j. If 
B >_ C^, purchase one unit of system j. 

If B < C^, select the largest 

A, 0) (t) 

— = such that C(i) < B. If no units can 

C (T ) 

be purchased, stop the procedure. 

Step 4 . If the item purchased in Step 3 is item j, 

next compute , . \ , . \ 

A 1 1 5 ( t ) Ag ( t ) 
max {max [ — g - p - p -, c (j) ) 

i^j 

The next unit is assigned to the i tme where the maximum is 
taken on (if the item can be purchased with the funds avail- 
able). Continue this procedure until adding an additional 
unit would exceed the budget available. 

Numerical Example . Suppose the budget available is 
B = $6,000, n = 4 systems, and the parameters are as given 
by Table 3. 

Table 4 shows the order in which the allocation procedure 
selects items for support. The procedure stops with $5,700 
spent since the remaining $300 cannot purchase any additional 
spares. The circled numbers indicate the order of selection. 
For example, the first item selected was item 2, then item 3, 
then i tern 1 , etc . 
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Table 3. Parameters for Resource Allocation Example. 



i 


1 


2 


3 


4 


C (i) 


$1 ,000 


$800 


$900 


$1 ,200 


MTBF 


30 days 


40 days 


50 days 


30 days 


MTTC 


5 days 


4 days 


2 days 


5 days 


P (i) 


0.8 


0.8 


0.9 


0.8 



Table 4. Cost-Benefit Ratios and Selection Sequence 



i 


1 


2 


3 


4 


A i ( t ) 

TT 


© 


© 


© 


© 


0.000348 


0.000457 


0.000404 


0.000290 


A 2 ( t ) 


© 


© 






TUT 


0.000207 


0.000210 


0.000164 


0.000173 


A,(t) 










cTiT 


0.000065 


0.000053 


0.000020 


0.000054 



The resource allocation procedure described above could 
easily be modified to reflect some sort of weighting scheme 
such as military essentiality codes if it is desired to 
consider the importance of the item in the allocation procedure. 
In that case, we would consider the ratios 

CUT 

where is the item weighting factor. The procedure would 

then work just as described above. Again, no claim of 

optimality can be made for this marginal allocation procedure, 

but practical experience with such schemes shows that they 
usually produce good results. 
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V. THE AVAILABILITY PROBLEM WITH CONSIDERATION 
OF SYSTEM CONFIGURATION 



In all of the preceding material we have assumed that 
the system consisted of a single component. This was done 
so that we could ignore the system configuration. In the 
real world problem, we would certainly not want to ignore 
system configuration, for the support levels of components 
of systems should depend on what happens to the system when 
the components fail. To determine this, one must examine 
the system configuration. The supply support procedure 
(COSAL) presently in use by the U.S. Naval Supply System 
totally ignores system configuration. In fact, it ignores 
the system itself, looking only at the components of those 
systems. Its contribution is that of piece part support. 

In this chapter we examine the impact that system configura- 
tion might have on operational availability and through an 
example we compare the allocations with the configuration 
considered and with the system configuration ignored. 

Complex weapon systems usually consist of several major 
subsystems or components which are actually the items supported 
aboard ship. We therefore consider each system as broken 
down into its "supportable" components. We now view our 
objective as one of maximizing "system" availability. We 
seek to achieve our objective by supporting the components of 
that system. We confine our discussion to weapon systems 
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consisting of components arranged in series, parallel or 
simple combinations of both. We assume that each weapon 
system is a coherent system of order n, and its subsystems 
are Type-1 items. The analysis for Type-2 items is identical. 

Let A s (t) and P S y S be the operational availability of 
the weapon system, at time t and its specified availability 
goal, respectively. Let A^^k.t) be the incremental system 
availability at time t obtained by adding one extra spare 
(from k-1 to k) of component i. 



Step 1 . Determine the system availability, A S y S (t), 
from the component availability equations assuming zero 
spares for every component. 

Step 2 . Compute max { — ^ — }. (The max is taken 



over all components in the system.) If the maximum is taken 
on for i=j and the budget is sufficient to purchase one unit 
then allocate C(j) dollars for component j. If not, select 
that component with the largest ratio that can be purchased 
with the available funds. 



Step 3 . Compute max (max [ — ^ q !{) — 



A (j) (2,t) 
i sys v * ’ 

— mr~ 



} and 



i 7* j 

continue as in the allocation procedure of the last chapter 
until no additional components can be purchased. 



Step 4 . The shipboard allowance list for spare stock 
consists of the spares determined by this procedure. 
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Numerical Example. Consider the missile system of Figure 6 



extracted from [7]. 



A 

A 



B C 

— 4—* 



i m 



i 



:i 



] — L 





U T ~ T” 

fj°n n s neC "i 

lSyitein_ j, ^Sy_stem J 



- 5 







1^* “T 

•Adjacent , 

^Systems 



Figure 6. Weapon System Configuration. 



Suppose that functionsA, B, and C are necessary for a 
successful intercept, and that al least one missile must be 
launched. The reliability block diagram for system function 
D (target intercept) is shown in Figure 7. 
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Figure 7. System Reliability Block Diagram. 



Consider that the weapon system has 11 independent 
subsystems whose parameters are given in Table 5. 



Table 5. Missile Subsystem Parameters 



System 

(1) 


MTBF ! 
(Days ) 


MTTC 
(Days ) 


Unit Price 
($> 


System 

(i) 


MTBF 
(Days ) 


MTTC 
(Days ) 


Unit Price 
($) 


1 


50 | 


i | 


j 5,000 


7 


30 


1 


10,000 


2 


60 | 


1 


7,000 


8 


30 


1 


10 ,000 


3 


100 1 

1 j 


4 


4,000 


9 


40 


2 


5,000 


4 


1 50 1 


2 


5,000 


1 0 


30 


1 


10,000 


5 


i 90 I 

I 


2 


2 ,000 


1 1 


30 


1 


10,000 


6 


100 


2 


5,000 











Using seri es -paral 1 e 1 reliability calculations, the opera- 
tional availability of the missile system is easily determined 
to be: 



A sys (t) ' A k' ) ( t >-''P ) (t)-A(3)( t) . A (4) (t) . A (5) (t) . {A (6) (t) . 

• £ 0 4 R fi 



[A k ?)(t) + A k S) (t) 3 + Ap ) (t).[A,[ 10) (t) + A^ 11} (t)]} 
7 o k 9 k 10 k ll 
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= 0.80 and a budget 



For a system availability goal of P 

s y b 

of B = $76,000 and the component parameters given in Table 5, 
Table 6 shows the detailed computational steps with the 
allocation sequence indicated by the circled values. The 
column headings are the incremental benef i t- to-cost ratios 

J- L. 

generated by adding an extra unit of the i component to 
the present listing of spares shown in column 1. The numbers 
in the vectors listed in column 1 denote the number of spares 
of each of the components. For example ( 2 , 1 ,0 , 1 ,0 ,0 ,0 , 1 ,0 ,0 ,0 ) 
means that there are presently 2 units of component 1, 1 unit 
of component 2, 1 unit of component 4 and 1 unit of component 
8. In the example, the funds were used up before the avail- 
ability goal was attained. The final spare parts determination 
for the missile system is ( 3 ,2 , 1 , 3 ,2 , 0 ,2 ,0 ,0 ,0 ,0 ) and the system 
availability at t = 45 is A s y s (45) = 0.7097. For comparison 
we used the same example and determined the spare parts listing 
using the allocation scheme described in the last section of 
the previous chapter. Table 7 gives the detailed computational 
results for that procedure. As before, the numbers in circles 
represent the sequence of allocations. The final spare parts 
listing when the system configuration is ignored is given by 
( 2 , 1 , 0 ,2 ,2 ,0 , 1 , 1 , 1 , 1 , 1 ) . With this allowance list the system 
availability is given by 

A S y S ( 45 ) = 0.5183 
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Table 6. Alocation Procedure for Missile System 
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Indifferent purchase among 7th, 8th, 10th and 11th systems. 



Table 7. Allocation Procedure for Missile System (Component Oriented). 
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We note that this availability is far below what was obtained 
when the system configuration was explicitly considered even 
though the latter procedure spent $1000 more. This is only 
one example out of many possibilities, but it does point out 
the importance of considering the system configuration in 
spare parts allocation schemes. We cannot claim that the 
improvement would always be as dramatic as that witnessed in 
this example but there must always be some improvement. 

After all, the objective is to maximize system availability, 
not the component availabilities. Therefore, it only makes 
sense to focus on the system and not the components. 
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VI. SUMMARY AND CONCLUSIONS 



In this chapter summary and conclusions are made 
concerning the discussions and results presented in this 
report. Following the conclusions, suggestions are given for 
extension and enrichment of this study. 

We have derived a mathematical expression to compute 
the operational availabilities for a standby spares system 
having no repair capability when there are finitely many 
spare parts for support. By employing a Markovian approach 
which, in turn, permits us to work with linear homogeneous 
differential equations with constant coefficients, we have 
determined a formula for calculating the operational avail- 
ability for a standby spares system with a repair capability. 
We have discussed decision rules for determining if a repair 
capability should be purchased, a decision rule for deter- 
mining the number of spares for the two systems, and finally 
resource allocation schemes for components within a system. 

We have shown that the resource allocation scheme should 
consider the system configuration when making decisions about 
those items that should be supported. Certainly, if the 
attainment of a specified level of system availability is 
the objective, then resources should be allocated to support 
those components which are most critical, not necessarily 
those with the highest rates of failure. In our example, we 
witnessed a potential money savings of about 45% considering 
system configuration as compared to the one which neglects 
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the system configuration. These results support the argument 
that the Supply System should become more maintenance and 
system oriented and less piece parts oriented. 

In describing the operational availability of a given 
system (non-mai ntai ned system) it is necessary to specify (1) 
the equipment failure process, (2) the replacement mechanism, 
(3) the system configuration, and (4) the state in which 
system is to be defined as failed. For maintained systems, 
we should also consider the repair mechanism. The simplest 
hypothesis from a mathematical viewpoint is to assume that 
the equipment (component) failure processes, replacement 
processes and repair processes have exponential distributions. 
These assumptions make it relatively easy to determine the 
operational availability of a system as a function of the 
availabilities of its components when there are limited 
numbersof spare parts for support. These operational avail- 
ability functions, in turn, allow us to compute the correct 
operational availability of a complex system by considering 
its system configuration. As demonstrated in Chapters IV and 
V, these functions can then be used as decision criteria for 
the determination of spare part listings for support, such as 
the COSAL. When the dollars available for constructing a 
shipboard parts list are not sufficient to achieve specified 
levels of operational availability for all systems, the 
resource allocation scheme should consider the marginal 
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increase in system operational availability per dollar invested 
in order to maximize the system operational availability. It 
is also emphasized that system configuration should be 
considered throughout the resource allocation process. With- 
out these considerations, the limited resources might be spent 
in an inefficient way or perhaps additional funds would be 
expended to achieve the same level of operational availability 
which could otherwise have been obtained less expensively. 

In modeling the operational availability functions, 
exponential distributions have been assumed not only for the 
failure process, but also for the replacement and repair 
processes. This might not be the case in many applications. 
Since, in many cases, the use of the exponential failure law 
can be justified, this might not be critical for the component 
lifetimes. However, for the replacement and repair mechanisms, 
the assumption is probably unrealistic. Relaxation of the 
assumption results in complicated expressions for operational 
availability. Even with the simplest model shown in (9), it 
is very laborious to compute A k (t) values, when k gets large. 
This suggests the importance of developing good approximation 
methods. Even though this paper tries to solve basic issues 
associated with the shortcomings of current COSAL procedure, 
many problems such as the definition of operational avail- 
ability, military essentiality, system configuration, 
availability goals, and mission scenarios still require much 
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work before many significant changes can be made to the COSAL 
procedure. It is recommended that further study along these 
lines be conducted. 

It is hoped that this study will be helpful to those 
people who endeavor to improve the current COSAL procedure 
of the Naval Supply Systems Command. 
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